Persistence of Manifolds in Nonequilibrium Critical Dynamics 
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We study the persistence P(t) of the magnetization of a d' dimensional manifold (i.e. the prob- 
ability that the manifold magnetization does not flip up to time t, starting from a random initial 
condition) in a d-dimensional spin system at its critical point. We show analytically that there are 
three distinct late time decay forms for P{t) : exponential, stretched exponential and power law, 
depending on a single parameter C, = {D — 2-\-r])/z where D = d — d' and r?, z are standard critical 
exponents. In particular, our theory predicts that the persistence of a line magnetization decays as 
a power law in the d = 2 Ising model at its critical point. For the d = 3 critical Ising model, the 
persistence of the plane magnetization decays as a power law, while that of a line magnetization 
decays as a stretched exponential. Numerical results are consistent with these analytical predictions. 

PACS numbers: 05.70.Ln, 05.50.+q, 05.70.Jk 



Following a rapid quench of a spin system from infi- 
nite temperature to zero temperature, domains of com- 
peting ground states form and grow with time [1]. In 
the standard nonconserved dynamics a given spin inside 
the sample flips only when a domain wall passes through 
it, which happens rather rarely. As a result, persistence, 
i.e. the probability that the spin remains unflipped up to 
time t decays slowly as a power law, P{t) ^ at late 
times [2]. In contrast, if the spin system is quenched to its 
critical temperature Tc, any given spin fluctuates rapidly 
due to the finite temperature and the persistence of a sin- 
gle spin has an exponential tail at late times. However, 
for a quench to Tc, the persistence of the total magneti- 
zation of the sample (as opposed to that of a single spin) 
again decays as a power law, P{t) ~ t"^', where the ex- 
ponent 9g has been argued to be a new nonequilibrium 
critical exponent [3]. The global persistence has since 
been studied in a wide variety of systems [4] . 

In a d-dimensional sample, a single spin is a 0- 
dimensional manifold. On the other hand, the global 
magnetization corresponds to summing over all the spins 
of the full sample which is a d-dimensional manifold. 
To interpolate between these two limits, it is natural 
to study the persistence of the magnetization of a d'- 
dimensional manifold with < d' < d. In the two lim- 
its d' = and d' = d, the persistence has very differ- 
ent asymptotic decay, respectively exponential and power 
law. The following question then naturally arises: if one 
tunes the manifold dimension d' from d' = to d' = d, 
how does the asymptotic behavior of persistence changes 
from an exponential (d' = 0) to a power law decay 
(d' = d)? Does this change occur abruptly at some inter- 
mediate value of d', or is there an intermediate regime of 
d' where the behavior is neither exponential nor power 
law, but something in between? In this Letter we address 
this interesting issue and we show analytically that in- 
deed there is a novel intermediate regime of d' where the 
persistence has a stretched exponential tail. Our main 



results can be summarized in terms of a single number 
C = {D — 2 + r])/z, where D = d — d' is the co-dimension 
of the manifold and ry and z are the standard critical ex- 
ponents: z is the dynamical exponent that describes the 
temporal growth of the correlation length ^(t) ^ t^^^ and 
T] describes the power law decay of the equilibrium spin- 
spin correlation function at Tc, {(j){0)(j){r)) ^ 
for large separation r. Depending on the value of ^, the 
persistence of the magnetization of a d'-dimensional man- 
ifold at Tc has the asymptotic behavior. 



Pit) 



^-e(d' .d) 

exp (— oii'') 
exp {-bit) , 



C<0 
< C < 1 

C>i, 



(1) 



where the exponent 0(d', d) depends on d' and d and ai, 
hi are constants. Strictly speaking, in the intermediate 
regime Q < Q < \ we will show that, exp (— a2t^lni) < 
P{t) < exp (— ai<^) for large t. We will derive the results 
in Eq. (1) within the mean-fleld theory (valid for d > 4), 
in the n —^ oo limit of the 0{n) model, followed by a 
general scaling theory. 

Two speciflc applications of our general results are as 
follows: (i) Consider the persistence of the line magne- 
tization (d' = 1) in the d = 2 Ising model at Tc. Us- 
ing d' = 1, d = 2, ?7 = 1/4 and z w 2.172, one gets 
^ = {d - d' - 2 + r])/z w -0.3453 < 0. Hence Eq. 
(1) predicts a power law decay for the persistence of the 
line magnetization; (ii) For the d = 3 Ising model, using 
1] w 0.032 and z w 2, one finds that while for the plane 
magnetization (d' = 2), ( « —0.484 < 0, for the line 
magnetization (d' = 1), < C ~ 0.016 < 1. Our Eq. 
(1) then predicts a power law decay for the persistence 
of the plane magnetization, but a stretched exponential 
decay for the line magnetization. Numerical simulations 
for small samples in 2 and 3 dimensions are consistent 
with these analytical predictions. 

Our starting point is the standard Langevin equation 
for the vector order parameter (j) = (0i, . . . , 0„), 
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dt(j)i = V'^<j)i - r(j)i - (u/n)0^ ^<P'3+Vi, (2) 

where ?7(x, t) is a Gaussian white noise with mean zero 
and correlator {r]i{x,t)r]j{x' ,t')) = 2SijS'''{yi—x')S{t—t'). 
The magnetization of a d'-dimensional manifold will 
be defined by the vector field, ipii^d'+i, ■ ■ ■ ^ Xd,t) = 
J d>i{x, t) HiLi obtained by integrating the order pa- 
rameter over the d! directions. Here L denotes the length 
of the sample along any direction. For a vector order pa- 
rameter, we define the persistence of the manifold to 
be the probability that any given component of the man- 
ifold magnetization, say Vi; does not change sign up to 
time t. Since all the components of the spin are equiva- 
lent, henceforth we will drop the subscript i of i/'i for con- 
venience. Note that the magnetization is a field over 
the remaining D = d — d! dimensional space whose co- 
ordinates {xd'+i-: • • • ) ^d) will be relabelled by the vector 
r = (n, r2, . . . ro) for convenience. The observation that 
allows us to make analytical predictions for the persis- 
tence of the magnetization 'ip{f, t) is that it is a Gaussian 
variable at all finite times. This follows simply from the 
fact that 'ip{f'^ t) is a sum of L'^ random variables which 
are correlated but only over a finite correlation length 
^{t) ~ t^^^. Thus in the thermodynamic limit when 
t^^^ ^ L, the central limit theorem asserts that tp{f, t) is 
a Gaussian field. Hence its persistence is determined by 
the autocorrelation function C{t\,t2) = {'tjj{f,ti)il){f,t2)) 

We start with the moan-field theory, valid for (i > 4, 
where we set u = 0, and also r = (at the critical 
point) in Eq. (2). Next wc integrate the Langcvin equa- 
tion over the d' space directions and then solve the re- 
sulting linear equation in the Fourier space. Defining 
ip{k, t) = J ip{r, t)e^'''^d^r, it is easy to compute the two- 
time correlation function, 

(V'(fc,ii)V'(-^,t2))=A(^)e-'='(*i+*^) 

+ J_(e-fe^|ti-*2|_e-fe^(ti+t.)j^ 

where A(fc) = 0)V'(— fc, 0)) is taken to be a con- 
stant (for uncorrelated initial condition). At late times, 
the initial condition dependent term becomes negligi- 
ble and it is sufficient to retain only the second term 
on the right-hand side of Eq. (3). The autocorrela- 
tion function is then obtained by integrating over the k 
space, C{ti,t2) = J d'^ke-''^''\i){k,ti)i}{-k,t2)) where 
we have introduced a soft ultraviolet cut-off a. One finds 

Citut2) = PA [(ti +t2 + a^f^ - -t2\+ a^f^] , 

(4) 

where 2(3 = {2 — D)/2 and A is an unimportant constant. 

Consider first the case D < 2, i.e., /3 > 0. In this 
case there is no need for the ultraviolet cut-off since 



{xp'^{r,t)) = C{t,t) does not diverge at any finite time 
even if one puts a = 0. Setting a = in Eq. (4) we 
find that the correlation function C{ti,t2) is still non- 
stationary. However, one can render it stationary by 
employing a well known trick [2], where one introduces 
a normalized Gaussian process, X = tp/^/^tp^ which, 
when observed in the logarithmic time T = Int be- 
comes a stationary Gaussian variable with the correlator, 
C(T) = (X(O)X(T)) = cosh(T/2)^^- |sinh(T/2)p/5. In- 
terestingly, the same Gaussian correlator also appears in 
the context of the persistence of the Edwards- Wilkinson 
type rough interfaces [6]. It is well known [5,2] that for 
such a stationary Gaussian correlator (decaying exponen- 
tially for large T) , the persistence of the process also de- 
cays exponentially for largo T, P{T) ^ exp(— ^T). In 
terms of the real time t = , this indicates a power law 
decay, P{t) ~ t~^, with persistence exponent 9. For this 
particular correlator, the exponent has been studied in 
great detail in the context of the interface problem [6] 
and the exponent is known to depend continuously on 
the roughness parameter /3 = (2 — £))/4 > 0. 

For D > 2, on the other hand, ono needs the ultravio- 
let cut-off a explicitly in order to keep {tp'^{f, t)) = C{t, t) 
finite. In that case, the appropriate scaling limit is 
ti,t2 oo, but keeping their difference \ti — ^2! fixed. In 
this limit, Eq. (4) reduces to a stationary correlator in the 
original time variable, C{ti,t2) ~ (|ti -^2! +a^)"^^~^^/^ 
that decays as a power law for large \ti — t2|. To cal- 
culate the persistence of Gaussian stationary processes 
with an algebraically decaying correlator is nontrivial. 
However, there exists a powerful theorem due to Newell 
and Rosenblatt [7], which states that if the station- 
ary correlator decays as C{t) ~ t~°' with a > for 
large time difference t = \ti — t2\, then the persis- 
tence P{t) (probability of no zero crossing between ti 
and t2) of such a process has the following asymptotic 
behaviors: (i) P{t) ^ cxp{—Kit) if a > 1 and (ii) 
exp(-i4:2t"lnt) < P{t) < exp(-if3r) if < a < 1, 
where the Ki's are constants. Applying this theorem to 
our problem, we find that P{t) ^ cxp(—Kit) for D > 4 
and exp(-if2t(^-')/'lnt) < P{t) < ex^i-K^t^'^-^^'^) 
for 2 < D < A. In the borderline case i3 = 4, there will 
be an additional logarithmic correction. 

Combining these results for D < 2 and D > 2 and not- 
ing that z = 2 and ?7 = within the mean-field theory, we 
find that the explicit exact results for the mean-field the- 
ory derived above are just the special cases of the general 
result in Eq. (1) provided one uses the mean-field value 
C = (£> - 2)/2 in Eq. (1). 

The mean-field theory is valid for d > 4. In order to 
access the physically relevant dimensions d < 4, we now 
consider another solvable limit, namely the n — > 00 limit 
where Eq. (2) becomes 

^t(j>i = V^(j>i-[r + S{t)](P^+^J^, (5) 
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where S{t) = u{(f>f) has to be determined self- 
consistently. The critical point corresponds to r + 
5(oo) = 0. This self-consistent determination of S{t) can 
be done using standard techniques [3] and one finds that 
at late times, S{t) S{oo) - (4 - d)/U for 2 < d < 4. 
Substituting this result into Eq. (5), summing over the d' 
directions, solving the resulting equation in the Fourier 
space and finally integrating over the k space (as in 
the mean-field theory), we finally arrive at the following 
autocorrelation function for the manifold magnetization 
ilj{f,t) in 2 < d < 4, 

(6) 

where ti <t2, Ai is a constant and a represents the soft 
ultraviolet cut-off as before. 

For I? < 2, as in the mean- field theory, one can set the 
cut-off a = and the resulting nonstationary correlator 
in Eq. (6) can be reduced to a stationary correlator for 
the normalized process X = il)/ ^y^tp^ in the logarithmic 
time T = Int, 



A{T) 



where fi = {d - 2)/2, 2/3 = (2 - D)/2, B[m,n] is the 
standard Beta function and B[m, n, x] = dyy"^^^{l — 
y)"'~^. Since the stationary correlator in Eq. (7) decays 
exponentially for large T, A{T) - exp[-(d + i:>- 2)r/4], 
one concludes [5,2] that the corresponding persistence 
also decays exponentially for large T, P{T) ~ exp(— 0T) 
and hence as a power law in the original t = variable, 
P{t) ~ t'^ . Determining the exponent analytically is 
still a challenging task. However one can make progress 
in the limit of small co-dimension D ^ Q. Note that for 
D = 0, i.e. for the global persistence, Eq. (7) becomes 
a pure exponential A{T) = cxp[-(d - 2)T/4] for all T, 
indicating that the process is Markovian [3]. One then 
finds P(T) ~ exp(--6'oT) with = {d - 2)/4 [3]. For D 
nonzero but small, one can expand the correlator in Eq. 
(7) around the Markov process [D = 0) and then use a 
perturbation theory result [8] to calculate 9 to first order 
in D. We get 9 = 0o + Dd^Id/TT + OiD^) for all 2 < rf < 4, 
where Id is a complicated d-dependent integral. For spe- 
cial values of d, this integral simplifies [9]. For example, 
for d = 4, we get, 6* = 1/2 + (2V2 - l)i:'/4 + 0{D'^) and 
for d= 3, = 1/4-1- 0.183615... Z)-t-0(D2). 

For I? > 2, it is evident from Eq. (6) that one needs to 
keep a nonzero cut-off a in order that the integral does 
not diverge at t\ = t2- In this case, for large ti, the 
dominant contribution to the integral comes from the 
regime — > ii. It is easy to see that in the limit when 
ti,t2 are both large with their difference \ti — t2\ fixed, 
the autocorrelator in Eq. (6) reduces to a stationary one. 



C{ti,t2) ~ - t2\ + a2)-(^-2)/2, where Bi is an 

unimportant constant. One can then invoke the Newell- 
Rosenblatt theorem [7] once again to conclude that for 
large t, the persistence P{t) ^ exp(— Kit) if D > 4 
and exp(-K2t(^-'^/'lnf) < P{t) < eM-Kst^'^'^^^^) 
for 2 < Z) < 4, where the Kj's are constants. Combining 
these results we thus find that the n ^ oo results for 
the persistence of manifold magnetization in 2 < d < 4 
are also compatible with our general results in Eq. (1) on 
noting that C = (-D — 2)/2 since r] = and z = 2 within 
the large n limit. 

Taking hints from the two solvable cases above we 
now construct a general scaling theory valid for all d > 
2. The two-point correlation function of the order pa- 
rameter, at the critical point, has the generic scaling 
form, {(j){0,ti)(j){x,t2)) - x~(-''~^+'^'>F{xt^^^"',t2/ti) for 
large distance x and large times ti, t2, where r] and z 
are the standard critical exponents defined in the in- 
troduction. This indicates that in the Fourier space, 
_(^(K, ii)^(-K, t2)) ~ K-(-^-'i^G{KtY%t2/ti), where K 
is a c?-dimensional vector conjugate to x. The manifold 
magnetization tl^ is obtained by summing the order pa- 
rameter (j) over d' directions. This is equivalent to putting 
Ki = along the i = 1, . . . ,d' directions. One then 
obtains the scaling behavior of the two point correla- 
tor of the manifold magnetization, {ip(k,ti)ip{—k,t2)) ~ 
k~^'^~^'> g(ktY^' ,t2/ti), where k is now a £> = d — rf' di- 
mensional vector. For example, within the mean- field 
theory, the scaling function g{x,y) = exp[— x^jl — y\] — 
cxp[— a:;2(l + y)], as evident from Eq. (3) after dropping 
the A dependent term at late times. The autocorrelation 
function C{ti,t2) ~ (■0(r, ii)V'(r, ^2)) is then obtained by 
integrating over k, 



C{tut2) 



d'^k 



g{ktY'',t2/ti)e~ 



a 



(8) 



where a is the soft ultraviolet cut-off as before. 

Consider first the case when D — 2 + -q < 0. One 
can then set the cut-off a = (since the integral 

in Eq. (8) is convergent at the upper limit), and one 



obtains C{ti,t2) ~ t 



-{D-2+-q)/z 



f{t2/ti) in the limit 



ti,t2 00 with t2/ti arbitrary. Note that the function 

f{x) ~ x-^-/^ for large x such that C{ti,t2) ~ i^^"^"" 
for t2 S> tl where Ac is the standard autocorrelation 
exponent [10,11,1]. This nonstationary Gaussian cor- 
relator can then be reduced, as before, to a station- 
ary one for the normalized variable X = -0/ (ip^) in 
the logarithmic time, T = Int and one gets, A{T) = 
(X(O)X(r)) = exp[(I? - 2 + 7?)r/2^]/(e^)//(l). Since 
A{T) ~ exp[-(Ac-(D-2-hr?)/2)T/z] for large T, it 
follows, as before, that the persistence P{T) ^ cxp{—6T) 
for large T. This means that the persistence decays as a 
power law in the original time variable t = , P{t) ~ 
for large t. 
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In the complementary case D — 2 + rj > 0, the inte- 
gral in Eq. (8) is, for ti = t2, divergent near the up- 
per limit without the cut-off. Hence one needs to keep 
a nonzero and then the appropriate scaling limit is ob- 
tained by taking ti, t2 both large keeping their difference 
1^1—^2! fixed but arbitrary. Then one can replace the scal- 
ing function g{kt-^ , t2/ti) in Eq. (8) by another function 
giik\ti — 12\^^^) of a single scaling variable, as in the two 
previous solvable cases. Performing the integral, one then 
finds C{h,t2) - |ti - ^21"^^"'+"^/" for \h - t2\ > a\ 
This correlator is stationary and decays as a power law. 
Invoking the Newell-Rosenblatt theorem once more, we 
find that P{t) decays exponentially for {D — 2 + ri)/z > 1 
and as a stretched exponential for 0< {D — 2 + ri)/z < 1. 
Combining this with the result for _D — 2 + ry < outlined 
in the previous paragraph, gives our general result in Eq. 
(1) on defining ( = {D - 2 + r])/ z. 
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FIG. 1. The three outer curves represent the persistence of 
the line magnetization in d = 2 critical Ising model as a func- 
tion of time for system sizes L = 63, 95 and 127. The two 
curves in the middle represent the persistence of the plane 
magnetization in d = 3 critical Ising model for system sizes 
L = 15 and 31. For the same system sizes, the two leftmost 
curves represent the persistence of the line magnetization in 
d = 3. 

To test the analytical predictions (i) and (ii), we have 
done preliminary Monte Carlo simulations for d — 2 
and c? = 3 Ising model for small size lattices with pe- 
riodic boundary conditions. The results are summarized 
in Fig. 1. The systems were evolved using heat-bath 
Monte Carlo dynamics at their bulk critical couplings: 
fsTc = [ln(l + \/2)]/2 in d = 2 and ~ 0.221656 in 
d = 3. For d = 2, the data for the persistence of the 
line magnetization for system sizes L = 63, 95 and 127 
(in each case the data was averaged over 1000 samples) 
shows a power law decay P{t) ~ t^^ for t <^ and 
crosses over to a faster decay for i ^ . A fit to the lin- 
ear part of the data on the log-log plot gives an estimate 



9 « 0.72. Similarly the persistence for the plane magne- 
tization in d = 3 shows a power law decay, P{t) ^ 
with 6 ~ 0.88, estimated from small lattice sizes L = 15 
and 31. The estimates of these exponents are only rough 
and may shift a little with larger lattices. The persis- 
tence of the line magnetization in d = 3, in contrast, 
has a much faster decay. Note that our theory predicts 
a decay, P{t) ~ exp(— oit^) where ( ~ 0.016. Such a 
small stretching exponent is difficult to determine from 
the small size lattices due to the strong finite-size effects 
at late times. All we can say is that the present data 
for the line magnetization in d = 3 is consistent with a 
decay of persistence that is faster than a power law but 
slower than an exponential. More extensive simulations 
and a somewhat sophisticated finite size scaling analysis 
are required to pin down the precise value of the stretch- 
ing exponent [9] . 

In summary, we have shown that the persistence of 
submanifolds in a critical system decays with time in 
a manner that depends on the dimensions d and d' of 
the system and manifold respectively. The crossover be- 
tween the power-law decay observed for the 'global' per- 
sistence (d' — d) and the exponential decay of the 'local' 
persistence (d' = 0), occurs via an intermediate regime 
of stretched exponential decay with stretched exponent 
( = [d-d' -2 + T])/z, for < C < 1- This latter behavior 
is predicted for the line magnetization of the d = 3 crit- 
ical Ising model. Numerical simulations are consistent 
with the analytical predictions. 
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